Eio, BEio and Arithmetical Chaos in Superstring Cosmology 
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It is shown that the never ending oscillatory behaviour of the generic solution, near a cosmological 
singularity, of the massless bosonic sector of superstring theory can be described as a billiard motion 
within a simplex in 9-dimensional hyperbolic space. The Coxeter group of reflections of this billiard 
is discrete and is the Weyl group of the hyperbolic Kac-Moody algebra Eio (for type II) or BEio 
(for type I or heterotic), which are both arithmetic. These results lead to a proof of the chaotic 
("Anosov") nature of the classical cosmological oscillations, and suggest a "chaotic quantum billiard" 
scenario of vacuum selection in string theory. 
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One crucial problem in string theory is the problem 
of vacuum selection. It is reasonable to believe that 
this problem can be solved only in the context of cos- 
mology, by studying the time evolution of generic, in- 
homogeneous (non-SUSY) string vacua. In this vein, 
it has been recently found |Q,0] that the general solu- 
tion near a spacelike singularity (t — > 0) of the massless 
bosonic sector of all superstring models (D = 10 IIA, 
IIB, I, HE, HO), as well as that of M-theory (D = 11 
SUGRA), exhibits a never ending oscillatory behaviour 
of the Belinskii-Khalatnikov-Lifshitz (BKL) type ||. In 
this letter, we analyze in more detail the asymptotic dy- 
namics (as t —* 0) of this oscillatory behaviour. We find 
that the evolution of the scale factors and the dilaton at 
each spatial point can be be viewed as a billiard motion 
in some simplices in hyperbolic space H 9 , which have 
remarkable connections with hyperbolic Kac-Moody al- 
gebras of rank 10. 

The central idea of the BKL approach is that the var- 
ious points in space approximately decouple as one ap- 
proaches a spacelike singularity (t — > 0). More precisely, 
the partial differential equations that control the time 
evolution of the fields can be replaced by ordinary dif- 
ferential equations with respect to time, with coefficients 
that are (relatively) slowly varying in space and time. 
The details of how this is done are explained in || for 
pure gravity, and in for the graviton-dilaton-p-form 
systems relevant to superstring/M-theory. 

We shall focus in this letter on the dynamical be- 
haviour of the metric and the string dilaton and recall 
first the relevant equations from ^|). To leading or- 
der, the metric (in either the Einstein frame or the string 

frame) reads g^ dx^ dx v = ~N 2 (dx°) 2 + J^Li a] M 2 , 
where d = D— 1 denotes the spatial dimension, and where 
lo 1 (x) = A (x) dxi is a d-bein whose time-dependence 
is neglected compared to that of the local scale factors 
a,. It is convenient to work with the 10 field variables 
fi — 1, ...,10, with, in the superstring (Einstein- 
frame) case, P l = — In en (i = 1, . . . , 9), and /3 10 = — (p 
where ip is the Einstein-frame dilaton. [In M-theory there 
is no dilaton but fj, = i = 1, . . . , 10. In the string frame, 
we define [3$ = — \n(y/gse~ 2 ®) and label fi = 0, . . . , 9.] 



We consider the evolution near a past (big-bang) or fu- 
ture (big-crunch) spacelike singularity located at t = 0, 
where t is the proper time from the singularity. In the 
gauge N — — Jg (where g is the determinant of the 
Einstein- frame spatial metric), i.e. in terms of the new 
time variable dr = —dt/^/g, the action (per unit comov- 
ing volume) describing the asymptotic dynamics of (3^ as 
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In addition, the time reparametrization invariance 
(i.e. the equation of motion of N in a general 
gauge) imposes the usual "zero-energy" constraint E — 
G^{d(3^ /dT){d(3 v /dr) + V(/?") = 0. The metric 
in field-space is a 10-dimensional metric of Lorentzian 

signature !-+•••+. Its explicit expression depends 

on the model and the choice of variables. In Af-theory, 

G% ddl, d^ M = Ei° =1 (d^ M ) 2 - (e;°=i d^ M )\ while in 

the string models, d(3% df3 v s = £i=i (d/3 l s ) 2 - {dj3%) 2 
in the string frame. Each exponential term, labelled by 
A, in the potential U(/3 M ), Eq. (|J), represents the ef- 
fect, on the evolution of {g^ v ,Lp), of either (i) the spa- 
tial curvature of c/y ("gravitational walls"), (ii) the en- 
ergy density of some electric-type components of some p- 
form ("electric p-form wall"), or (iii) the energy 

density of some magnetic- type components of A^,,,^ 
("magnetic p-form wall"). The coefficients Ca are all 
found to be positive, so that all the exponential walls in 
Eq. (^) are repulsive. The CU's vary in space and time, 
but we neglect their variation compared to the asymp- 
totic effect of wa(P) discussed below. Each exponent 
— 2wa{0) appearing in Eq. (^) is a linear form in the 
/3 M : u>a (0) — wah • The complete list of "wall forms" 
WAiP), was given in [jjj for each string model. The num- 
ber of walls is enormous, typically of the order of 700. 

At this stage, one sees that the r-time dynamics of 
the variables /3 M is described by a Toda-like system in a 
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Lorentzian space, with a zero-energy constraint. But it 
seems daunting to have to deal with ~ 700 exponential 
walls! However, the problem can be greatly simplified 
because many of the walls turn out to be asymptotically 
irrelevant. To see this, it is useful to project the motion of 
the variables onto the 9-dimensional hyperbolic space 
H 9 (with curvature —1). This can be done because the 
motion of (3^ is always time-like, so that, starting (in 
our units) from the origin, it will remain within the 10- 
dimensional Lorentzian light cone of G^ v . This follows 
from the energy constraint and the positivity of V . With 
our definitions, the evolution occurs in the future light- 
cone. The projection to H 9 is performed by decomposing 
the motion of (3^ into its radial and angular parts (see 
and the generalization ||; see also 0] for recent 
comments on the covariance of the chaos obtained in the 
billiard approximation). One writes /3 M = +pj^ with 
p 2 = -G^u 0*j3 v ,p>O and G^ Y 7" = -1 (so that 7" 
runs over H 9 , realized as the future, unit hyperboloid) 
and one introduces a new evolution parameter: dT = 
kdr/p 2 . The action (m becomes 



S = k dT 
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where d*y 2 — G^ dj" is the metric on H 9 , and where 
V T = k- 2 p 2 V = k- 2 C AP 2 exp(-2puM( 7 )). 

When t — ► + , i.e. p — > +00, the transformed potential 
Vr(p, 7) becomes sharper and sharper and reduces in the 
limit to a set of p-independent impenetrable walls located 
at 10^4(7) = on the unit hyperboloid (i.e. Vr = when 
Wa{j) > 0, and Vr = +00 when 104(7) < 0). In this 
limit, d hip/dT becomes constant, and one can choose 
the constant k so that d\np/dT — 1. The (approxi- 
mately) linear motion of (3^ (r) between two "collisions" 
with the original multi-exponential potential V{(3^) is 
thereby mapped onto a geodesic motion of f(T) on H 9 , 
interrupted by specular collisions on sharp hyperplanar 
walls. This motion has unit velocity (dj/dT) 2 — 1 be- 
cause of the energy constraint. In terms of the original 
variables the motion is confined to the convex domain 
(a cone in a 10-dimensional Minkowski space) defined by 
the intersection of the positive sides of all the wall hy- 
perplanes wa(P) = and of the interior of the future 
light-cone ft 1 * j3 u = 0. 

A further, useful simplification is obtained by quo- 
tienting the dynamics of (3^ by the natural permuta- 
tion symmetries inherent in the problem, which corre- 
spond to "large diffcomorphisms" exchanging the vari- 
ous proper directions of expansion and the corresponding 
scale factors. The natural configuration space is therefore 
M. d /Sd, which can be parametrized by the ordered multi- 
plets (3 1 < 1 < ■ ■ ■ < f3 d . This quotienting is standard 
in most investigations of the BKL oscillations || and can 
be implemented in K d by introducing further sharp walls 
located at f3 % = (3 l+1 . Note that the natural permutation 
symmetry group is different in M-theory (where it is S10), 
and in the D = 10 string models (Sg), and would be still 



smaller in the successive dimensional reductions of these 
theories. However, there is a natural consistency in quoti- 
enting each model by its natural permutation symmetry. 
Indeed, one finds that, upon dimensional reduction, there 
arise new (exponential) walls, which replace the missing 
permutation symmetries in lower dimensions || . Finally 
the dynamics of the models is equivalent, at each spa- 
tial point, to a hyperbolic billiard problem. The specific 
shape of this model-dependent billiard is determined by 
the original walls and the permutation walls. Only the 
"innermost" walls (those which are not "hidden" behind 
others) are relevant. 

We have determined the set of innermost walls for all 
string models. The analysis is straightforward || and we 
report here only the final results, which are remarkably 
simple. Instead of the O(700) original walls we find, in 
all cases, that there are only 10 relevant walls. In fact, 
the seven string theories M, IIA, IIB, I, HO, HE and the 
closed bosonic string in D = 10 P], split into three sep- 
arate blocks of theories, corresponding to three distinct 
billiards. The first block (with 2 SUSY's in D = 10) is 
B2 = {M, IIA, IIB} and its ten walls are (in the natural 
variables of M-theory (3^ 



B 2 :wf ] (/3) = -?+f? +l (i = l,...,9) 
4o(/3) =/3 1 +/3 2 +/3 3 - 
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The second block is B\ = {I, HO, HE} and its ten walls 
read (when written in terms of the string-frame variables 
of the heterotic theory a 1 = (3g, a° = 0s) 



61 :w^ ] (a) = a 1 , wf ] (a) = -a'" 1 + a l (i = 2, . . . ,9), 



w w ( a ) = a° - a 7 - a 8 ~ a 9 . 



(5) 



The third block is simply Bq = {D = 10 closed bosonic} 
and its ten walls read (in string variables) 



Bo : w [ " ] (a) = a 1 + a 2 , w l " l (a) = -a^ 1 + a\i = 2,..., 9), 



wf ] (a) = a - a 7 - a 8 - a g 
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In all cases, these walls define a simplex of H 9 which is 
non-compact but of finite volume, and which has remark- 
able symmetry properties. 

The most economical way to describe the geometry of 
the simplices is through their Coxeter diagrams. This 
diagram encodes the angles between the faces and is ob- 
tained by computing the Gram matrix of the scalar prod- 
ucts between the unit normals to the faces, say r[™' = 

-w^ where Wi = Wi / ■ Wi, i = 1, . . . , 10 labels the 
forms defining the (hyperplanar) faces of a simplex, and 
the dot denotes the scalar product (between co-vectors) 
induced by the metric G M „ : Wi ■ Wj = Wi^ Wj V for 
Wi{0) = Wi fl f3 fJ '. This Gram matrix does not depend 
on the normalization of the forms Wi but actually, all 
the wall forms Wi listed above are normalized in a nat- 
ural way, i.e. have a natural length. This is clear for 
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the forms which are directly associated with dynamical 
walls in D = 10 or 11, but this can also be extended 
to all the permutation-symmetry walls because they ap- 
pear as dynamical walls after dimensional reduction [|S|. 
When the wall forms are normalized accordingly (i.e. 
such that y i d y namical ^ exp(— 2 Wi(j3)), they all have a 

squared length wf^ • wf^ = 2, except for ■ vdp = 1 
in the B\ block. We can then compute the "Cartan ma- 
trix" , a'™' = 2 wf^ ■ I wf^ ■ iwj n ' , and the corresponding 
Dynkin diagram. One finds the diagrams given in Fig. 
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FIG. 1. Dynkin diagrams defined (for each n = 2,1,0) by 
the ten wall forms (/3 M ), i = 1, . . . , 10 that determine the 
billiard dynamics, near a cosmological singularity, of the three 
blocks of theories B 2 = {M, IIA, IIB}, Bi = {I, HO, HE} and 
Bo = {D = 10 closed bosonic}. The node labels 1, . . . , 10 
correspond to the form label i used in the text. 

The corresponding Coxeter diagrams are obtained 
from the Dynkin diagrams by forgetting about the norms 
of the wall forms, i.e., by deleting the arrow in BEio. As 
can be seen from the figure, the Dynkin diagrams asso- 
ciated with the billiards turn out to be the Dynkin di- 
agrams of the following rank-10 hyperbolic Kac-Moody 
algebras (see [fTo)): Ei , BEi and DEi (for £> 2 , B\ and 
Bq, respectively). It is remarkable that the three bil- 
liards exhaust the only three possible simplex Coxeter 
diagrams on H 9 with discrete associated Coxeter group 
(and this is the highest dimension where such simplices 
exist) H]. The analysis suggests to identify the 10 wall 
forms wf 1 ^^), i = 1, . . . , 10 of the billiards B2, B\ and 
Bq with a basis of simple roots of the hyperbolic Kac- 
Moody algebras E10, BE10 and DE10, while the 10 dy- 
namical variables /3 M ,// = 1,...,10, can be considered as 
parametrizing a generic vector in the Cartan subalgebra 
of these algebras. It was conjectured some time ago JlJ] 
that E10 should be, in some sense, the symmetry group 
of SUGRAn reduced to one dimension (and that DE10 
be that of type I SUGRA10, which has the same bosonic 
spectrum as the bosonic string). Our results, which in- 
deed concern the one-dimensional reduction, a la BKL, 
of M/string theories exhibit a clear sense in which E10 
lies behind the one-dimensional evolution of the block B2 
of theories: their asymptotic cosmological evolution as 



t — > is a billiard motion, and the group of reflections 
in the walls of this billiard is nothing else than the Weyl 
group of E10 (i.e. the group of reflections in the hyper- 
planes corresponding to the roots of E10, which can be 
generated by the 10 simple roots of its Dynkin diagram). 
It is intriguing - and, to our knowledge, unanticipated 
(see, however, |13 )- that the cosmological evolution of 
the second block of theories, B\ = {I, HO, HE}, be de- 
scribed by another remarkable billiard, whose group of 
reflections is the Weyl group of BEio- The root lattices 
of E10 and BE10 exhaust the only two possible unimodu- 
lar even and odd Lorentzian 10-dimensional lattices Jio| . 

A first consequence of the exceptional properties of the 
billiards concerns the nature of the cosmological oscilla- 
tory behaviour. They lead to a direct technical proof 
that these oscillations, for all three blocks, are chaotic in 
a mathematically well-defined sense. This is done by re- 
formulating, in a standard manner, the billiard dynamics 
as an equivalent collision- free geodesic motion on a hyper- 
bolic, finite-volume manifold (without boundary) A4 ob- 
tained by quotienting H 9 by an appropriate torsion-free 
discrete group. These geodesic motions define the "most 
chaotic" type of dynamical systems. They are Anosov 
flows E4J , which imply, in particular, that they are "mix- 
ing" . In principle, one could (at least numerically) com- 
pute their largest, positive Lyapunov exponent, say X^ n \ 
and their (positive) Kolmogorov-Sinai entropy, say h^. 
As we work on a manifold with curvature normalized to 
— 1, and walls given in terms of equations containing only 
numbers of order unity, these quantities will also be of 
order unity. Furthermore, the two Coxeter groups of E10 
and BE10 are the only reflective arithmetic groups in H 9 
p"T[ so that the chaotic motion in the fundamental sim- 
plices of E10 and BE10 will be of the exceptional "arith- 
metical" type |lq ]. We therefore expect that the quan- 
tum motion on these two billiards, and in particular the 
spectrum of the Laplacian operator, exhibits exceptional 
features (Poisson statistics of level-spacing,. . .), linked to 
the existence of a Hecke algebra of mutually commuting, 
conserved operators. Another (related) remarkable fea- 
ture of the billiard motions for all these blocks is their 
link, pointed out above, with Toda systems. This fact 
is probably quite significant, both classically and quan- 
tum mechanically, because Toda systems whose walls are 
given in terms of the simple roots of a Lie algebra enjoy 
remarkable properties. We leave to future work a study 
of our Toda systems which involve infinite-dimensional 
hyperbolic Lie algebras. 

The present investigation a priori concerned only the 
"low-energy" (E <§; (a')" 1 / 2 ), classical cosmological be- 
haviour of string theories. In fact, if (when going to- 
ward the singularity) one starts at some "initial" time 
to ~ (dp/dt)^ 1 and insists on limiting the application 
of our results to time scales \t\ (a') 1 / 2 = t s , the to- 
tal number of "oscillations" , i.e. the number of colli- 
sions on the walls of our billiard will be finite, and will 
not be very large. The results above show that the 
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number of collisions between to and t — > is of order 
7V coU ~ lnr ~ Ln(ln(i /f)). This is only iV coU - 5 if 
t corresponds to the present Hubble scale and t to the 
string scale t s . However, the strongly mixing properties 
of geodesic motion on hyperbolic spaces make it large 
enough for churning up the fabric of spacetime and trans- 
forming any, non particularly homogeneous at time to, 
patch of space into a turbulent foam at t = t s . Indeed, 
the mere fact that the walls associated with the spatial 
curvature and the form fields repeatedly rise up (during 
the collisions) to the same level as the "time" curvature 
terms ~ t~ 2 , means that the spatial inhomogeneities at 
t ~ t s will also be of order tj 2 , corresponding to a string 
scale foam. 

Our results on the 62 theories probably involve a 
deep (and not a priori evident) connection with those 
of Ref. [H on the structure of the moduli space of M- 
theory compactified on the ten torus T 10 , with vanishing 
3-form potential. In both cases the Weyl group of E 10 
appears. In our case it is (partly) dynamically realized as 
reflections in the walls of a billiard, while in Ref. (l6) it is 
kinematically realized as a symmetry group of the moduli 
space of compactifications preserving the maximal num- 
ber of supersymmetries. In particular, the crucial E-type 
node of the Dynkin diagram of E10 (Fig. [I]) comes, in our 
study and in the case of M-theory, from the wall form 

[2] 



{(3) = (3\j + + 0m associated with the electric 



energy of the 3-form. By contrast, in |16| the 3-form 



is set to zero, and the reflection in w\q comes from the 
2/5 duality transformation (which is a double T dual- 
ity in type II theories), which exchanges (in M-theory) 
the 2-brane and the 5-brane. As we emphasized above, 
dimensional reduction transforms kincmatical (permuta- 
tion) walls into dynamical ones. This suggests that there 
is no difference of nature between our walls, and that, 
viewed from a higher standpoint (12-dimcnsion ?), they 
would all look kinematical, as they are in |16|. By anal- 
ogy, our findings for the B\ theories suggest that the Weyl 
group of BE10 is a symmetry group of the moduli space 
of T 9 compactifications of {I, HO, HE}. 

Perhaps the most interesting aspect of the above anal- 
ysis is to provide hints for a scenario of vacuum selec- 
tion in string cosmology. If we heuristically extend our 
(classical, low-energy, tree-level) results to the quantum, 
stringy (t ~ t s ) and/or strongly coupled (g s ~ 1) regime, 
we are led to conjecture that the initial state of the uni- 
verse is equivalent to the quantum motion in a certain 
finite volume chaotic billiard. This billiard is (as in a 
hall of mirrors game) the fundamental polytope of a dis- 
crete symmetry group which contains, as subgroups, the 
Weyl groups of both Ei and BE10 |p7| . We are here as- 
suming that there is (for finite spatial volume universes) a 
non-zero transition amplitude between the moduli spaces 
of the two blocks of superstring "theories" (viewed as 
"states" of an underlying theory). If we had a descrip- 
tion of the resulting combined moduli space (orbifoldcd 
by its discrete symmetry group) we might even consider 



as most probable initial state of the universe the funda- 
mental mode of the combined billiard, though this does 
not seem crucial for vacuum selection purposes. This 
picture is a generalization of the picture of Ref. [[[9) and, 
like the latter, might solve the problem of cosmological 
vacuum selection in allowing the initial state to have a 
finite probability of exploring the subregions of moduli 
space which have a chance of inflating and evolving into 
our present universe. 
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